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Abstract 

We consider a set of coupled underdamped ac-driven dynamical units exposed to 
a heat bath. The coupling scheme defines the absence/presence of certain symme- 
tries, which in turn cause a nonzero/zero value of a mean dc-output. We discuss 
dynamical mechanisms of a dc-current appearance and identify current reversals 
with synchronization/desynchronization transitions in the collective ratchet's dy- 
namics. 
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1 Introduction 

The ratchet effect, i. e. the possibility to obtain the directed transport by 
using zero-mean perturbations only, has induced notable scrutiny over the 
last decade [1]. Initially, most studies have been focused on noisy overdamped 
models have been inspired by a molecular motors realm [1]. Then the ratchet's 
approach has been applied to a broad class of physical systems in which inertia 
effects are essential [2]. The examples are Josepshon junctions [3], cold atoms 
systems [4], and mechanical engines [5]. 

Recently it has been shown that the ratchet idea can be viewed as a part 
of a general symmetry-breaking approach [6]. This approach is based on the 
analysis of all relevant symmetries which have to be broken in order to fulfill 
necessary conditions for a dc-output appearance. The formalization of sym- 
metry analysis for one particle's dynamics has been addressed in Ref.[7]. 
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In the present paper we aim at collective rectification effects which arise in a 
set of coupled single ratchet units. Although various examples of interacting 
ratchet systems have been proposed already in the context of molecular motors 
[8], here the emphasis is put on the weak- noise underdamped limit. We show 
that a coupling scheme determines a set of certain symmetries for the ratchet's 
collective and hence defines necessary conditions for a dc-current generation. 
Dynamical mechanisms of a current rectification are connected with a coher- 
ence between units, which depends not only on a coupling scheme, but also 
on a strength of interactions. 



2 Coupling schemes and symmetries 

Let us consider a set consisting of N identical dynamical units, x = {xi,i = 
1,...,N}, that are linearly and symmetrically coupled. The coupling scheme 
is described by some graph which can be encoded in the symmetrical N x N 
binary matrix Gij, Gij = Gji. The equations of motion are the following: 

N 

Xi = -aii + F(xi - x°,t - i°) + c^2GijH(xi,Xj,Xi,Xj) + (1) 

i=i 

where F(x + L, t) = F(x,t + T) = F(x, t) is the double periodic force function, 
H is linear over all arguments coupling function, and c is the strength of inter- 
actions. The stochastic terms are mutually independent delta-correlated 
Gaussian white noises, = crSij8(t — s), where a is the noise inten- 

sity. We assume also that the force function, F, is the same for all the units, 
but the phases, xf and tf, can be different for different units. Finally, we are 
interested in the mean dc-current, 

1 N 

J = lim —YxAt). (2) 

t-oo Nt-, y ' 

Following the symmetry analysis ideology [6-7], in order to determine neces- 
sary conditions for a dc-output appearance, we have to check whether there 
exist symmetry transformations which allow to generate out for each trajec- 
tory of the system (1) another one with a reversal velocity. The presence of 
the white noise, does not change the symmetry properties of the system. 
Moreover, the coupling to a heat bath leads to an effective exploration of the 
whole phase space and produce an averaging in a case of several coexisting 
attractors [6,9]. 

For the one-particle case, N — 1, a transformation of interest has to involve 
a change of the sign of x (and thus change of the current direction). It allows 



2 



also some shifts in time and space domains [6-7]: 

S single ■ x -> -x + A, t -> t + r, i G R. (3) 

In the case of several coupled ratchets a symmetry operation should be per- 
formed in the global coordinate space H N and can also involve a permutation 
between different units. The corresponding symmetry operation can be de- 
scribed as the linear transformation, 

Snetwork ■ X — > -Sx + A, t — > t + T, XG^, (4) 



where SV,- is & N x N binary matrix with only one non-zero element at each 
row and line, and A = = 1, ...,N is the vector of shifts. The matrix S 

encodes permutation between units. 

The properties of the permutation matrix SV,- (and the existence of such a 
matrix at all) are strongly depend on the structure of connections among 
units. Below, using simple examples, we will illustrate this statement. 

Breaking symmetries by connections. Let us first consider the case of two par- 
ticles (rotators) in the standing-wave potential with the modulated amplitude 
[10], coupled by the linear spring: 



xi = -ax x + F(xi,t) + c(x 2 - xi) + £i(£) (5) 
x 2 = -ax 2 + F(x 2 — x ,t — t ) + c(x l - x 2 ) + (6) 

where F(x,t) = sin(x) sin(u;t). 

In the uncoupled case, c = 0, both the systems posses symmetries of the type 
(3): 



Si : x -> -x + ir,t -* t + T/2 (7) 
S 2 :x^-x + 7T + x ,t^t + T/2 + t , (8) 

with Xq = t = for the first rotator, Eq.(5). The symmetry transformations, 
Eqs.(7-8), are independent for each rotators. The whole system, Eqs.(5-6), is 
symmetric with respect to the transformation Si x S 2 for any choice of to an d 
xq. So, the mean dc-output for the uncoupled case is zero (line(l) in Fig.la). 

In the case of coupled particles, c > 0, abovementioned transformations, , 
Eqs.(7-8), should been conjugated and independent spatial and temporal shifts 
are forbidden now. For x ^ kn and to 7^ T/2, the connection breaks both the 
symmetries and we can expect on nonzero current appearance (see also Ref. 
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[11] for another example of an overdamped dimer with an additive driving 
force). For the set of parameters a = 0.1, w = 0.3, Xq = tt/2, a = 0.01 
and to — T/4, the symmetry violation is realized by the asymmetrical limit 
cycle with the negative winding number (line(2) in Fig. la). So, in this case 
the connection between the units destroy all the symmetries and leads to the 
dc-current generation. 




-3 3 

X 

Fig. 1. (a) The dependence of x\(t) versus t for different coupling schemes (see text 
for details) for the parameter values a = 0.1, w = 0.3, xo = tt/2 and to = T/A 
and noise intensity a = 0.01. The dependencies for another rotators from the set 
are the same due to the homogeneity of the system, Eqs.(5-6) and Eqs.(9-10); (b) 
Poincare sections for the first rotator for the third variant of the coupling scheme. 
The coupling to the heat bath leads to the averaging over two symmetry-related 
chaotic attractors (white dots) with opposite mean velocities. 

Restoring symmetries by connections. Let us now consider the dimer identical 
to the previous one, Eq.(5-6), but spatially shifted by ir (half of the period): 

x 3 = -ax 3 + F(x 3 + 7r, t) + c(x 4 - x 3 ) + £ 3 (t) (9) 
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X4 = -ax^ + F(x 2 - x + 7r, t - t ) + c(rc 3 - x A ) + £4(2)- 



(10) 



Both the systems, Eqs.(5-6) and Eqs.(9-10), can be transformed from one 
to another by the simple coordinate shift. Thus, the system in Eqs.(9-10) 
produces the same mean current, as the system in Eqs.(6-7). 

As the next step we couple both the systems by additional connections ((1) <-> 
(4), (2) <-> (3)) (see the inset in Fig. la). The generalized symmetry transfor- 
mation, Eq.(4), can be identified with the following permutation matrix: 

^1 0^ 



S. 



1 
10 
10 



(11) 



and the shift vector A = {n}. 

Here, despite to the previous case, the introduction of additional connections 
leads to the appearance of the new symmetry and we may expect on the nul- 
lification of the dc-output (line(3) in Fig. la). For the above set of parameters 
we found that the symmetry is realized in the phase space by two symmetry- 
coupled chaotic attractors with opposite mean velocities (see Fig. lb). It is easy 
to check, that any other type of connection with two additional links does not 
restore the system's symmetry. 



3 Synchronization and current reversals 



While the presence or the absence of the dc-output is clearly connected to 
the absence/presence of the symmetry transformation, the dc-current value 
is determined by dynamical mechanisms. From the previous studies [12] it 
is known that global properties of a collective dynamics are defined by a 
coherence between units. Thus, we can expect that an efficiency of the current 
rectification is closely related to a degree of synchronization within a ratchet's 
collective. 

As a natural example we consider the model of iV globally coupled array of 
underdamped Josepshon junctions (JJ), subjected to an ac-current E(t) [13]. 
The equation for the superconducting phase difference Xi across the single 
junction is 



Xi = -axi + sin(xi) + E(t) + c(xi - (x) A ) + 



'12) 
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Fig. 2. The dependence of the mean dc-output, J, and the degree of decoherence, 
D, versus the strength of interaction c for the system of N = 50 coupled units from 
Eq.(12). The values of parameters are a = 0.1, w = 1, a = (circles) and a = 0.03 
(triangles). 

where (x)a — jjJ2iLi x i{t) is the instantaneous mean array current. As the 
driving force we used the two-harmonics combination, E(t) = cos{ut) + 
cos(ut + 7r/2), which ensures that all the relevant symmetries are broken [7]. 
The numerically obtained dependence of the array dc-output on the strength 
of interaction c is shown in Fig. 2(a). For the set of parameters a = 0.1, uj = 
1,0" = 0, the dependence demonstrates the presence of two successive current 
reversals, at c ~ 0.03 and c ~ 0.33. In order to understand dynamical mech- 
anisms of these events, we introduced the mean decoherence, which we define 
as 



where (■■■)a,t means the averaging over the array and over the one period 
of ac-driving. The first result from the comparison of both the dependencies 
(Fig.2(b))is that the current reversal at c ~ 0.33 is connected with the tran- 
sition to the regime of the complete synchronization, D = 0. 

Lets now track briefly a relation between the current reversal and the synchro- 
nization. Initially, in the uncoupled limit, c = 0, all JJ's work independently, 
each one as the single one-dimensional rectificator [1,6]. The asymmetry of 
the ac- force E(t) is realized through the chaotic attractor with the positive 
mean current (see Fig. 3a) [14]. The local phase space of each units, (xi,Xi,t), 
has identical structure, but the relative phases of different units are randomly 
distributed (depending on the initial conditions). 

After the introduction of a nonzero coupling, c > 0, some coherence between 
the units occurs. From a point of view of the dynamics of the single junction, 
this leads to changes of an attractor structure in the phase space R 3 (see 




(13) 
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Fig. 3. Poincare sections for a single junction from the array of N = 50 coupled JJ, 
Eq.(12), for (a) c = and (b) c = 0.27 (zero- noise case). Parameter values are the 
same as in Fig.2. 



Fig. 3b). This causes changes in a projection of the attractor's invariant density 
on the velocity subspace x,{ and, as a result, leads to the current reversal at 
c pa 0.03. The further increase of the interaction strength up to c ss 0.33 results 
in the complete synchronization and to a shrinking of the global attractor in 
R 2Ar+1 to the hyperplane R 3 ,( t). The global attractor has now 

dimension equal to 3. The attractor of the single unit now is the same as in the 
uncoupled limit but all the units have the same phase. The mean dc-output 
returns to its value at the limit c = 0. 

Here some analogy with a current reversal in an one-dimensional deterministic 
ratchet [14] can be drawn. In the one-dimensional case current reversals have 
been identified with tangent bifurcations from chaotic to regular (limit cycles) 
attractors [14]. In our case the current reversal at c ~ 0.33 corresponds to 
the transition "hyperchaos - chaos" [15], connected with the shrinking of a 
the system attractor in a global phase space R 27V+1 . It is interesting to note, 
that from a point of view of the local attractor of the single junction, this 
transition corresponds to a crisis, connected with expanding of the attractor 
in three-dimensional subspace (xi,Xi,t) (compare Fig. 3a and Fig.3b). 

The presence of a weak noise (which is equivalent to a weak coupling with a 
heat bath) can strongly suppress correlations between units and may lead to a 
delayed synchronization transition [16]. Due to a strong conjugation between 
the synchronization and the current rectification process, we can manage a 
noise-induced current reversal for a fixed strength of interaction (Fig. 2(b)). 
Thus, the presence of a thermostat allows to control the dc-output by changing 
a temperature of the system. 



4 Concluding remarks 

Finally, we have presented the symmetry approach to the problem of the 
collective current rectification by a set of coupled dynamical units. This idea 
can be used for a more general problem such as an obtaining of a non-zero value 
of some relevant mean ensemble characteristic, {A(x)) t . This characteristic 
can correspond, for example, to a mean magnetization of a spin lattice with 
a complex geometry [17]. The proposed collective ratchet's ideology may also 
be relevant in a context of a cooperative dynamics of neural networks, as an 
approach to a visual processing with directional selectivity [18]. 

Relations between symmetries of a collective ratchet and its coupling scheme 
on the one side, and relations between a topology of interactions and synchro- 
nization properties [19], on the other one, may open an interesting perspective. 
It has been found that in the thermodynamic limit, N — > oo, the synchro- 
nization is impossible for nearest-neighbor coupled dynamical networks if the 
number of sites connected to a given site, N c , is a finite fraction of the total 
number of sites N, = const [20]. On the other hand, in a case of small- world 
networks, i.e. sets with long (in a sense of a topological distance) connections, 
the synchronization can be achieved in the thermodynamic limit through a 
small fraction of distant connections [20]. This nontrivial effect provides a 
tool for a control of the dc-output in a massive collective ratchet by changing 
a small number of relevant connections only. 

This work has been supported by the Emmy Noether-Programm of the DFG 
under contract LU1382/1-1 
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